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INTRODUCTION 


Space missions proposed for the future call for large, flexible systems 
launched by the Shuttle (refs. 1 to 4). These systems can either operate at 
Shuttle altitudes or be moved to geosynchronous orbits or beyond. One means 
of moving these large systems without exceeding their low acceleration re- 
quirements is by the use of electric propulsion thrusters (ref. 5). Elec- 
tric propulsion is also proposed for scientific missions such as the comet 
explorers (refs. 6 and 7). 

These spacecraft will require large, multikilowatt power sources for 
thruster operations. The conventional design of electric propulsion systems 
call for power to be generated at about 200 volts and use power processing 
to convert this into voltages required for thruster operations (ref. 8). 
However, a study for the proposed Halley Comet mission concluded that the 
net payload could be increased if a "direct drive" electric propulsion sys- 
tem is used (ref. 9). In this direct drive system, the high voltages re- 
quired for thrusters ('^1.2 kV) are generated directly on the solar array. 
Similar improvements in payload capability should also be realized for other 
direct drive electric propulsion missions. 

It has been known since the early 70' s that the geomagnetic substorm 
environment can charge spacecraft surface causing electronic switching anom- 
alies and even satellite power system failures (e.g. , see refs. 10 to 12). 
This spacecraft-charging phenomenon is under study in a joint AF/NASA in- 
vestigation (ref. 13). The objective of this investigation is to develop 
techniques to control differential charging of geosynchronous satellites. 
Now, with the advent of much larger spacecraft that may have high-voltage 
surfaces exposed, these spacecraft-charging control techniques will be of 
greater importance to system designers. There have been several initial 
analytical treatments of large space system interactions with space en- 
vironments (refs. 14 to 17) and preliminary treatments of electric thruster 
interactions with spacecraft and environments (refs. 18 and 19). 

In this paper, the NASA Charging Analyzer Program (NASCAP) computer 
code (refs. 20 and 22) is used to compute voltage distributions around a 
Solar Electric Propulsion (SEP) spacecraft as it encounters an idealized 
geomagnetic substorm environment. Both a standard operating voltage and 
direct-drive voltage configuration are considered. The computations are 
presented first without thruster operations and then with a simplified, sim- 
ulated thruster-on representation for direct-drive configuration only. It 
should be stressed that these computations looking for possible areas of 
concern in the spacecraft design are exploratory. Future improvements in 
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modelling resulting from such initial studies will enhance the value of com- 
puter codes as design tools. 


ANALYTICAL MODEL 


NASCAP Description 

The NASA Charging Analyzer Program (NASCAP) hi s een described pre- 
viously in the literature (refs. 20 and 22) and is w>niy briefly summarized 
here. NASCAP is a quasistatic computational code, that is, it assumes that 
currents are functions of environmental parameters, electrostatic potential, 
and magnetostatic fields. It is capable of analyzing the charging of 
three-dimensional, complex bodies as a function of time for given space en- 
vironmental conditions. It includes consideration of dielectric material 
properties (e.g. , secondary emission, backscatter, photoemission, and bulk 
and surface conduction) and computes currents involving these materials in 
determining the surface voltages and potential distributions around the body 

In NASCAP, the body must be defined in terms of rectangular par- 
allelipipeds, sections of parallelipides or flat plates within a 17x17x33 
point grid. Seven separate conductors can be specified with the first con- 
ductor capable of floating with respect to space while the others can be 
biased with respect to the first. The environment can be defined in terms 
of single or double Maxwellian distributions (ref. 23) by specifying elec- 
tron and proton temperatures (in electron volts) and plasma number den- 
sities. The code outputs a variety of graphic displays showing the model 
used, the voltage distributions for given environments at specified times, 
and particle trajectories (if desired). 

The code can simulate low-density particle emitters either electrons or 
protons. Initial voltages, current density, and angles can be specified and 
the emission of either species computed for the field distribution around 
the body. If the charged particles escape, they are counted in the overall 
current balance and the potential of the floating conductor adjusted ac- 
cordingly. The code can not yet treat particle-to-particle interactions. 
Hence, an accurate treatment of ion thruster fluxes is not possible. 

Discharges on the surfaces of the body can also be simulated by NASCAP 
code. At a specified voltage difference between an insulator surface and 
the conductor beneath, a specified fraction of the charge deposited on the 
insulator will be transferred to the conductor beneath, the potential dis- 
tributions readjusted for the charge redistribution, and the charging proc- 
ess continued. This simple discharge simulation is a representation of bulk 
breakdown through an insulator and can be used, in design studies, to pin- 
point areas where large electrical stress can build up. 

Spacecraft Description 

The NASCAP model of the Solar Electric Propulsion spacecraft used in 
this report is shown in figure 1. It is representative of a 25-kilowatt 
version of such spacecraft. Each square in the NASCAP model is lOmxlOm to 
allow the 30mxl0m spacecraft to fit within the NASCAP grid point limitations 

The solar array wings are each 10mxl2m. They are modelled as thin flat 
plates with 0.013 cm (6 mil) silica cover glass on the cells and a O.OlO cm 
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(4 mil) Kapton aubstrate. The aluminum aquares on the arraya repreaenta the 
expoaed metallic interconnecta. Since NASCAP can not aimulate email gapa, 
theee interconnecta muat be lumped together aa ahown. Thia lumped area rep- 
reaenta 5X of the total aolar cell area and ia a reaaonable approximation to 
propoaed deaigna. Each wing ia divided into three areaa which are at apeci- 
fied voltagea relative to the apacecraft body. For the standard con- 
figuration, the three areaa are ^50 volta each (^150 volts total) and for 
the direct drive configuration the three areas are at +400 volts (+1200 
volta total). These voltage configurations are shown in figure 2. 

The spacecraft body is modelled in two parts: a simulated experiment 

deck and thruster system. The experiment deck is modelled as an octagon 4m 
across x Im thick. It is covered by 0.010 cm (4 mil) Kapton blanket with a 
2mx2m aquar> of exposed metal (instrument platform). The thruster system 
consists of a 2mx3mx4m section covered by 0.020 cm (b mil) optical solar 
reflectors (OSR) , 0.020 cm (4 mil) Kapton and aluminum (to simulate thruster 
areas), and a housekeeping section, 2mx2mxlm, covered by 0.010 cm (4 mil) 
Kapton. The experiment deck is connected to the thruster system by a me- 
tallic boom; both are assumed to be at spacecraft ground potential. The 
solar arrays are connected to the spacecraft body by fairly large metallic 
booms which are also assumed to be at spacecraft ground potential. 

In this model the criteria for occurrence of discharges is based on an 
edge voltage gradient breakdown concept (ref. 24) with the breakdown volt- 
ages determined from laboratory test data. For the Kapton surfaces break- 
downs are assumed to occur when electric fields exceed 1.3x10^ volts/cm, 
and for both the OSR and silica cover slides when the field exceeds 3x10^ 
volts/cm (ref. 23). The charge lost in a discharge is assumed to be 30)! of 
the charge scored on the insulator surface just prior to the discharge. 

This is based on a study that indicated Chat not all of Che scored charge in 
large area samples is removed in a discharge (ref. 26). 

Space Environment 

For this study an assumed geomagenCic substorm environment is used 
(fig. 3). This environment, utilizing single Maxwellian electron and proton 
temperatures, illustrates the transient nature of substorm environments and 
allows evaluation of spacecraft designs over a range of space conditions. 

The plasma densities chosen are large but reasonable. Since it has been 
shown Chat the charging rate is a strong function of Che number density 
(ref. 27), then use of the values in this model should promote rapid charg- 
ing of the surfaces and minimize computer time to expose Che areas of con- 
cern. 

Solar illumination is used in this study with Che angle of incidence at 
27° to the solar array normal. This simulates January or July conditions 
in orbit and enhances possible sun-sha<^.“ differential charging. 

DISCUSSION OF RESULTS 


Standard Configuration 


The response of selected surfaces on Che spacecraft model to Che sub- 
storm environment is shown in figure 4. In figure 4(a) the response of sur- 
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faces on Che experiment deck is shown. Since Che sunlit Kapton surfaces on 
Che deck generally follow Che spacecraft ground potentials, they are not 
plotted. In figure 4(b) Che response of surfaces on Che thruster system 
section are shown. The response of Che -ISO volt section of the solar array 
is shown in figure 4(c). The response of Che other sections of Che array 
are similar Co Chat of this section. 

The data in these figures indicates that the spacecraft ground responds 
to Che subsCorm by charging at a slower rate than the shaded insulators 
(e.g., see fig. 4(a)). This is to be expected since the capacitance of Che 
spacecraft as a whole is much larger Chan that of isolated insulators. This 
does indicate the development of differential charging on spacecraft Chat 
was assumed to be Che cause of spacecraft charging effects (ref. 28). Note 
Chat Che silica cover slides on the solar array are positive with respect Co 
Che interconnects (see fig. 4(c)). This is due Co Che photoand secondary 
emission characteristics of Che silica. 

Major discharges start between 7 and 8 minutes after encounter with the 
subsCorm. This is shown by the strong negative peaks in Che spacecraft 
ground potentials in all three figures (the solar array potential is fixed 
relative Co spacecraft ground and remains at this set bias throughout Che 
simulation). The negative peak results from Che NASCAP discharge simulation 
which requires that charge be transferred Co Che conductor beneath. The 
conductor voltage is then recalculated based on this new charge content. 

This is not a perfect simuilation of discharges since it is known Chat 
charge can be lost to space in Che discharge process (ref. 26). However, it 
does provide a graphic display of arcing rates. It can easily be seen Chat 
major arcs occur every 6 minutes in the first phase of the subsCorm (Tg * 

8 keV) and continues at a rate of once per 8 minutes in the second phase of 
the subsCorm (Tg <■ 6 keV). It should be noted Chat the silica cover glass 
on the solar array does not discharge: Che negative peaks shown simply fol- 

lows Che spacecraft ground (fig. 4(c)). 

The charging/discharging characteristics in both the spacecraft ground 
and insulator voltage traces appear to be very regular in the first phase of 
the subsCorm. However, in Che second phase irregularities begin to appear 
(note Che 30 to 33 min data of fig. 4(a)). These are due to breakdowns on 
only parts of the spacecraft. The voltage distributions around C)'e 
spacecraft at times of 31, 32, and 33 minutes are shown in figure 5. These 
figures show an end view through Che center of the spacecraft. These ir- 
regularities could be an effect due to the large size of the spacecraft 
since it did not appear in models of smaller satellites studied (ref. 27). 

It should be investigated further. 

In the third phase of the substorm (Tg > 3 keV) all discharges stop 
and the spacecraft comes into equilibrium with a spacecraft ground potential 
of about -1200 volts relative Co space. The insulator surfaces also stabi- 
lize at their respective values with the largest voltage difference being in 
the shaded Kapton at about 800 volts (fig. 4(a)). This substantiates that 
there is an electron temperature threshold in order for disch rges to oc- 
cur. This study and others (e.g., ref. 27) indicates Chat this threshold is 
at an electron temperature of about 5 keV. When Che spacecraft moves into a 
quiescent environment, it rapidly, but not drastically, looses its charge. 
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Dirtct-Drive Configuration 

Without Thrmter Operation ; In thia part of the study the spacecraft 
is assumed to move through the substorm with high voltages in the solar ar- 
ray and thrusters off. This is done to establish the behavior for com- 
parison with the standard configuration and with the simuilated thruster-on 
performance to be discussed in the next section. 

The charging behavior of the spacecraft in this configuration is shown 
in figure 6. In figure 6(a) the response of the experiment deck is shown. 
The response of the thruster system section and shaded OSR's is shown in 
figure 6(b). In figure 6(c), the response of the silica solar cell covers, 
interconnects and Kapton substrate of one of the 1200-volt sections are 
shown. 

The first observation is that spacecraft ground is driven about 800 
volts more negative than in the standard configuration. This is believed to 
be due to the high positive voltages in the array. The silica cover glass 
behavior is also different: it is now negative with respect to the inter- 

connects (figs. 4(c) and 6(c)). 

Discharges occur here as they do with the standard configuration but 
the regularity patterns have changed. In the standard configuration, large 
numbers of surfaces discharged at the same time setting up a regular 
charge-discharge sequence. In the direct-drive configuration the discharges 
start at about 10 minutes after substonn encounter (about 2 min later than 
in the standard configuration). Thereafter, it appears that breakdowns oc- 
cur piecemeal with a few surfaces discharging at each computational time 
step rather than massive numbers of surfaces discharging periodically. This 
indicates that, while the high voltages could cause more discharges, they 
could be less intense. 

The voltage distributions around the spacecraft for 16, 17, and 18 min- 
ute time steps are shown in figure 7. The principal difference that occurs 
is in the voltage distributions around the solar array. It is apparent that 
the high voltages in the array sections causes cor:.entrations whereas the 
standard configuration has fairly uniform fields at the solar arrays (com- 
pare figs. 5 and 7). 

When the third phase of the substorm (Tg ■ 3 keV) is encountered, the 
arcing again stops. The potentials appear to stabilize with the spacecraft 
ground at about -1800 volts and the highest voltage on the solar array wing 
at about -600 volts (fig. 6(c)). The insulators come to equilibrium values 
as determined by substorm conditions and leakage currents. When the sat~ 
ellite moves out of the substorm, the surface voltages all quickly relax 
(without discharges) to their normal, lower values. 

With Simulated Thruster Operations : As stated previously, the NASCAP 

code can not yet simulate ion thruster operations in a self-consistent man- 
ner. For this study the available proton and electron emitters are operated 
alternately for short periods of time to simulate the effects of plasma 
emission on spacecraft behavior. This simulation neglects 
part icle-to-particle interactions. However, this method should give a rea- 
sonable approximation to the spacecraft behavior if the alternating emission 
times are short compared to the charging times of the spacecraft surfaces. 

The following sequence is used in this study (see fig. 8). First, the 
spacecraft is allowed to charge negatively in the geomagnetic substurm for 
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lO cycles (570 sec). Then, while still in the first phase of the substons 
(Tg ■ 8 keV) , protons with a 1 mA bean current and an energy spread of 20 
eV to 2 keV are emitted for 3 milliseconds. This is to evaluate possible 
effects of accelerated and charge exchange plasma particles (ref. 29) on a 
negatively biased spacecraft. Then, electrons are emitted at a fixed beam 
current of 10 mA and energy of 1 keV. These values were selected to guar- 
antee that the electrons would leave the spacecraft and thus influence 
ground potentials. After 10 seconds, the electron emission is terminated 
and proton emission initiated for a total of 0.013 second. The electron 
emission is again started with the beam current rediced lu 0. 5 mA for the 
next 20 seconds. As the final phase of this simulation sequence, it was 
decided to evaluate the effects of constant accelerated electron emission 
only on spacecraft behavior. This is iccomplished by allowing the electron 
emitter to operate at 0. 5 mA beam cu;.-i.ent and 1 keV accelerating potential 
for the last 770 seconds. 

The results of this simulation are shown in figures 9 to 12. The com- 
putations of surface voltages for substorm conditions done are similar to 
the previous results. So, the discussion of results begins with the ini- 
tiation of proton emission at 570 seconds into the simulation. In figure 9 
the short operations with the protons being emitted are shown with an ex- 
panded time scale so that the effects could be illustrated. When the space- 
craft is negatively biased as it is just prior to emitter operations, proton 
emission first causes the ground potential to become less negative and then 
slowly became more negative. This is believed to be due to protons return- 
ing to spacecraft surfaces generating secondaries which escape. This 
changes the overall spacecraft current balance causing the shift in ground 
potential. The change in potential in turn changes the attractive potential 
for protons which modifies the current balance again causing the trend to 
more negative potential. It nnist be stressed that the overall system cur- 
rent balance is important and that this balance changes according to envi- 
ronment and emitter characteristics which influences the vehicle potential. 
Note that these effects are caused by relatively low currents - of the 1 mA 
of emitted protons only 50 to 80% return. 

When electrons are emitted, the spacecraft ground potential rapidly 
rises towards zero. It is interesting to note that insulators (such as the 
shaded Kapton experiment deck surface) follow spacecraft ground potential 
changes as predicted from ground-based simulation studies (ref. 30). Even- 
tually, the fixed beam current causes the spacecraft ground to go positive. 
At this point in time, electron emission is stopped and proton emission ini- 
tiated to relax the spacecraft ground potential towards zero again. The 
electron source is again initiated for 20 seconds and the spacecraft poten- 
tial remains close to zero volts. Finally, the proton only emission tends 
to drive the spacecraft slightly negative again. 

The NASCAP simiulations illustrate the effects of particle emissions on 
voltage distributions around the spacecraft (see fig. 10). Initially, the 
substorm has charged the spacecraft negative and the usual strong voltage 
distributions exist (see fig. 7). The initial proton emission does not 
change these distributions significantly (see fig. 10(a)). When electrons 
are emitted, spacecraft ground potential starts towards zero volts and the 
voltage distributions are significantly altered (figs. 10(b) and (c)). A 
voltage barrier (which reaches -500 volts relative to the spacecraft) builds 
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up. Had tha emitted electron energies been less than the barrier voltage, 
then emission would have been cut-off. This cut-off of low energy electron 
emission has been demonstrated in the ATS-S active charge control experi- 
ments (ref. 31). Proton emission tends to collapse this voltage barrier 
while the electrons tend to build it up. By alternating particle emissions 
in rather short time periods, the voltage distributions tend to remain 
rather constant (see figs. 10(d), (e), and (f)). This behavior is in agree- 
ment with ATS-6 active charge control experiments in which a plasma neutra- 
liser system was used to demonstrate that ground potentials could be con- 
trolled by emission of both ions and electrons (ref. 31). 

It should be noted that there are still strong electric fields that can 
exist during particle emissions (fig. 10). These are at the high voltage 
areas of solar arrays, shaded areas on the experiment deck, and shaded areas 
on the thruster system section. Discharges did occur in these areas during 
this simulation. However, they are not the major discharges that occur si- 
multaneously in a large number of cells and cause noticeable negative peaks 
in the spacecraft potential. 

Proton particle trajectories are interesting to view and a selection of 
these trajectories that occurred after spacecraft potentials were controlled 
is shown in figure 11. Under these conditions, there is no measurable cur- 
rent returning to the spacecraft surfaces - it all escapes to space. One of 
the supositions made to explain the ATS-6 active charge control experiment 
results was that the ions returned to the spacecraft neutralizing the sur- 
face charge which allowed the electrons to escape controlling the space- 
craft potential. Yet, in this case, the surfaces are negative but the pro- 
tons do not return. It could be that the positive charges effectively re- 
duce the barrier electric fields built up around the differentially charged 
spacecraft and this allows the electrons to escape. This is another topic 
for additional study. 

Since it has been found that alternating electron and proton emission 
appears to exercise beneficial control over spacecraft behavior in a charg- 
ing environment, then ion thrusters should be a practical means of providing 
this emission. Not only would they produce thrust but they could supply the 
relatively low currents of both electrons and ions required to control 
spacecraft potentials. 

Finally, the evaluation of high energy electron only emission is shown 
in figures 12(a) and (b). In these figures the time scale starts at 470 
seconds when the spacecraft is in the substorm and goes through the whole 
particle emission sequence. This discussion is limited to phase 7 simula- 
tion only. Electron emission with an accelerated beam does appear to con- 
trol spacecraft ground potentials. However, it does not seem to influence 
the insulator potentials. The results indicate that differential voltages 
build up to a point where large areas discharge. This occurs on the solar 
array and the experiment deck. 


CONCLUDING REMARKS 

The NASCAP computer code has been used to conduct preliminary computa- 
tions of the voltage distributions around large spacecraft in geomagnetic 
substorms. These spacecraft are possible configurations for future solar 
electric propulsion missions. Both a standard operating voltage configure- 
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tuion (♦ISO volt array) and a diracf-driva configuration (♦1200 volt array) 
are con'aidarad. The coaiputationa are aada with the thruatar*off for both 
configurations and with a thruater-on sinulation (eleetrona and protons 
mitted) for the direct*-drive configuration only. 

As with other satellites in geosynchronous orbit, it has been found 
that geonagnetic substoms can charge these spacecraft surfaces to strong 
negative potentials. Arc discharges in the standard and direct-drive con- 
figurations without operating eaitters arc highly probable. It appears that 
discharges in the direct-drive configuration are more infrequent, but less 
intense. 

It has been found that control of spacecraft voltages can be maintained 
by emission of relatively low currents; for example, 1 to 10 mA. It has 
also been found that alternating proton and electron emissions seems to con- 
trol voltage destributions around the spacecraft and minimise the volti^e 
barriers developed due to differential charging. This allows the spacecraft 
ground potentials to remain close to plasma potential. When the spacecraft 
ground is close to zero volts, protons all tend to escape. Since solar 
electric propulsion spacecraft have the inherent capability to emit both 
ions and electrons, they have a built-in active charge control technique. 

When only high energy electrons are emitted, the spacecraft grou^ po- 
tential can be controlled. However, differential charging does cause the 
development of a 500-volt barrier at the emitter. The insulator surface 
voltages are not significantly affected by this type of emission and even- 
tuallay large area discharges are observed. 

This study is a preliminary one, looking principally at possible ef- 
fects of the environment on large spacecraft. The simulation of ion thrust- 
ers used in this study is crude. It must be significantly improved to in- 
clude sitmiiltaneous, interactive emission before more positive conclusions 
can be reached on the effects of ion thruster operations on spacecraft sur- 
faces. 
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Figure 6. - Direct drive configuration behavior. 
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Figure?. - VoHage distriiHition around spacecraft. Direct-drive configuration; 
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Figure 8. - Simulated thruster operational sequence. 
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Figure 9. - Response o( spacecraft during simulated thruster operation. Direct-drive 
configuration in substorm. 
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Figure 10. - Voltage distribution around spacecraft Direct drive configuration in substorm: thruster 
simulation: equipotential lines in 100 V steps. 
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Figure 11. - Particle trajectories - protons. Beam current; 1 mA; energy range; 
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Figure 12. - Response to electron emission only. Direct-drive configuration in substorm. 
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